The X(3872) and other possible XYZ molecular states 
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Abstract. We perform a coupled channel calculation of the DD* and cc sectors in the framework of a constituent quark model. 
The interaction for the DD* states is obtained using the Resonant Group Method (RGM) and the underlying quark interaction 
model. The coupling with the two quark system is performed using the ^Pq model. The X(3872) is found as a molecular 
state with a sizable cc component. A comparison with Belle and BaBar data has been done, finding a good agreement. Other 
possible molecular molecular states are discussed. 
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INTRODUCTION 

In the last years a number of exciting discoveries of new hadron states, the so called XYZ mesons, have challenged 
our description of the hadron spectroscopy. Among them, one of the most mysterious states is the well established 
X(3872). It was first discovered by the Belle Collaboration in the J/xj/KK invariant mass spectrum of the decay 
K^n'^n^J /y/ [1]. Its existence was soon confirmed by BaBar [2], CDF [3] and DO [4] Collaborations. The world 
average mass is Mx = 3871 .2 ± O.SMeV and its width Fx < 2.3 MeV. The measurements of the X(3872) -> yj/\j/ 
decay [5, 6] implies an even C-parity. Moreover angular correlation between final state particles in the X(3872) 
Tt^Tt^J/yj/ decay measured by Belle [5] suggests that the 7^'' = 0++ and 7^'' = 0+^ may be ruled out and strongly 
favors the j'^^ = 1++ quantum numbers although the 2++ combination cannot be excluded. A later analysis by 
CDF Collaboration [7] of the same decay is compatible with the Belle results and concludes from the dipion mass 
spectrum that the most likely quantum numbers should be j'''- = 1++ but cannot totally exclude the J^^ — 2^+ 
combination . These conclusions were confirmed by a new CDF analysis of the decay X (3872) — ;r+7r^7/i^ followed 
by y/y/ — > 11^ jX^ excluding all the other possible quantum numbers at 99.7 % confidence level [8]. However the small 
phase space available for the decay X(3872) — > D^&^nP observed by Belle [9] discards the 7 = 2 leaving the 1++ 
assignment as the most probable option. 

In the 1++ sector the only well estabUshed state in the PDG [10] is the Xci[lP) with a mass M — 3510.66 ± 
0.07 MeV. The first excitation is expected around 3950Mey. In this energy region Belle has reported the observation 
of three resonant structures denoted by X(3940), 7(3940) and Z(3930). The last one was observed by Belle in the 
77— > DD reaction [11] and is already included in the PDG as the Xc2^P)- The X(3940) has been seen as a peak in 
the recoihng mass spectrum of J produced in e+e" collision. Its main decay channel is DD* [12]. The 7(3940) 
appears as a threshold enhancement in the J / \f/(0 invariant mass distribution of the B J/\j/coK decay [13]. Finally, 
CDF Collaboration has reported a new structure in the Z?+ — > J/ \j/<j)K^ decay at a mass of 4143 ±2.9 ± 1 .2MeV [14]. 
This state, called 7(4140) has some similarities to the 7(3940) state as far as decay channel is concerned. 

The decays measured for the X(3872) outlines a puzzling structure. In one hand the decays into nnJ/y/ and 
nnnJ / \lf through p and co mesons respectively suggests a sizable isospin breaking, incompatible with a cc structure. 
In the other hand the radiative decays into yJ/\j/ and yy/' suggest a sizable cc component. 



COUPLED CHANNEL CALCULATION 

In this work we present a coupled channel calculation of the 1++ sector including both cc and DD* states. The 
calculation is done in the framework of the constituent quark model of Ref. [15] widely used in hadronic spectroscopy. 



We start assuming a wave function given by 

a j8 

where \ \jfa) are cc eigenstates of the two body Hamiltonian, are cn (en) eigenstates describing the D (D) mesons, 
^Mai^) is the two meson state with j3 quantum numbers coupled to total j''^ quantum numbers and Xp (P) is the 
relative wave function between the two mesons in the molecule. The eigenstates of the C-parity operator are given by 
DD* =DD*±DD*. 

We use a phenomenological ^Pq model [16] to couple the two and four quark systems using the operator 



T = - 3 ^2/ £ y d^pd^p' 5^^\p + p') 



_ ,\ -1 C=l,/=0,S=l,/=0 

^i[^)bl{p)dl{p') 



(2) 



where H (v = jJ.) are the quark (antiquark) quantum numbers and Y = and 7 = ^ are dimensionless 

constants that gives the strength of the qq pair creation from the vacuum. The value of gamma is fitted to the 
Va(3770) -J> DD decay width. 

Using the wave-function from Eq. (1) with the two body eigenfunctions for all the mesons we apply the Resonant 
Group Method to obtain the dynamics in the two meson sector. We finnally end up with the coupled channel equation 

L / (n'p^^p' ,p)+yfp{P ,p))xp{p)p'dP=Exp\p') (3) 

where we include the and ^D\ DD* partial waves and 

y,t(p',p) = y:^>^^™i^ (4) 

P'p\ ^ ) Z- E-Ma 
is an effective interaction between the two mesons due to the coupling with intermediate cc states with 

{MM2p\T\Wa)=PVpo^{P)5^^\Pcm). (5) 



The cc probabiUties are given by 



^lllVap{P)Xp{P)P^dP (6) 



with the normalization condition 1 = La ka P + L/3 (Z/3 \Xp)- 

The results are given in Table 1 . Part A corresponds to the isospin symmetric case while part B shows the effect of 
isospin breaking in phase space. In both cases no bound state is found without couphng to the cc sector, neither in the 
7 = nor in the 7=1. This coupling generates a new state with an energy close to the DD* threshold. 

Having in mind that the ^Pq model is probably too naive and we might be overestimating the value of 7, we vary 
this parameter to get the experimental binding energy. The probabilities are given in part C of Table 1 . 



COMPARISON WITH THE DATA 

In order to compare the predictions of our model with the experimental data we use a Flatte-like parametrization of 
the DD* amplitude following Ref. [17, 18]. The differential cross section to final DD* states is given by 

dBrjB^KD^D*^) _ 1 I>o^ 

dE ~ |Z)(£)|2 ^'^ 

where ^ gives the branching to — )^ /rx(3872), r£,o£)*o(£) is the width calculated from the ^7}) model and 



D{E) = £-£/+^(r^o*o+ro+o-+r(£)) + ^(Ve/A'). 



(8) 



TABLE 1. Masses and channel probabilities for the three 
states in three different calculations. The first three states are 
found when we perform and isospin symmetric calculation with 
a value of y fit to the decay i|/(3770) DD. The second three 
states shows the effect of isospin breaking in the DD* masses. 
The last three states correspond to a value of 7= 0.19 that fits 
the experimental mass of the X(3872). The probability is shown 
as zero when it is less than 0.5 %. 
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FIGURE 1. Diagrams included in the quark rearrangement process DD* pJ/V- 



where r(£') accounts for the width due to other processes different from the opening of the near DD* threshold. 

The analysis of the B — > KX{3^12) Kn'^n^J/y/ data is more involved because we have to calculate the 
DD* n'^n~J I \\t transition amplitude. 

This can consistently be done in our formalism assuming that the process takes place through the DD* components 



of the X(3872) which decays into p//i/ and then into the final tt^ 
given by 



n J /\lf state. The decay width of the process is 



E 

JL 



dk- 







(9) 



(Mx-£p-£y/^)2 + -4 

The amplitude -^x^py/^ is calculated in our model by the rearrangement diagrams of Fig. 1, averaged with the 
DD* component of the X(3872) wave function. The rearrangement diagrams are calculated following Ref. [19]. The 
amplitude is given by 

i=a,d;j=b,h 

where 



(10) 
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FIGURE 2. Number of events for the decay B -> KD°D°k° measured by Belle (a) and for the decay B KD°D*° measured by 
BaBar (b). The solid and dashed lines shows the results from our model with and without the resolution functions as explained in 
the text. 
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FIGURE 3. Number of events for the decay B Kn^Tt .l/y/ measured by Belle (a) and by BaBar (b). The soUd and dashed 
lines shows the results from our model with and without the resolution functions as explained in the text. 



The spin-flavor-color matrix elements are taken from Ref. [19]. 

Once the decay width Fjj+jj-^/^ is calculated, the differential rate is given by 

dBr{B^Kn+n-J/\i/) _ 1 T^+^-j/^jE) 

dE " ^2^ |D(£)|2 • ^^^^ 

In order to compare with the experimental data we determine the number of events distributions from the differential 
cross section following Ref. [18]. In all reactions a background is taken into account modelled as in Ref. [18]. For 
the B KD^^ &^7t^^ the DD*° signal interferes with the background and so a phase (j)^""" = O" and (j)^"^"'' = 324^ 
have been introduced. Also the experimental branching ratio B{D*'^ — )• D'^nP) = 0.62 is introduced. We use a value for 
B§ = 3.5 10^"^ which is in the order of the one used in Ref. [18]. 

In Fig. 2 we compare our results with the B KD^D^n^ data from Belle (a) and B ^- KD^D*^ data from BaBar 
(b). The same comparison is done in Fig. 3 for the B — )• Kn^n~J/ Y data from Belle (a) and BaBar (b). In all figures 



the dashed lines shows the results without resolution functions. The solid line gives the result using the resolution 
functions as in Ref. [18]. All the resolution functions are those given by Belle [20] and BaBar [21] collaboration with 
the exception of the BaBar D^D*^ resolution where we use the prescription from Ref. [18]. 

We find a good description of the Belle B — )• KD^D^n^ data whereas the agreement is poor in the case of the BaBar 
data. It is important to notice that in the Belle analysis the mass of the X appears as 3E72MeV while in the BaBar data 
the resonance is located 3MeV above. The BaBar mass value does not coincide with the mass of the X obtained in our 
calculation which may be the reason for the disagreement. 

The B Kn^n^J /\)f data are equally well described for the Belle and BaBar experiments. In this case both 
Collaborations give similar values for the mass of the resonance, namely 3871.4Mey, which are in much better 
agreement with our result. 

Concerning other XYZ states,the X(3940), decaying to DD*, with mass M = 3942 ±9MeV is a good candidate to 
our state with a 88 % 1++ cc component and mass M — 3942Mey (see Table 1 part C). With respect to the other two 
states, it has been suggested that both, the 7(3940) decaying to J/\ifCO and 7(4140) decaying to J /\lf(p, are D*D* and 
D*D* hadronic molecules with 7^*- = 0++ or 2++ respectively. We have explore these channels but we have not found 
in principle any molecular bound state. 

As a summary, we have shown that the Z (3872) emerges in a constituent quark model calculation as a dynamically 
generated mixed state of a DD* molecule and Xci i?-P)- This structure allows to understand simultaneously the isospin 
violation showed by the experimental data and the radiative decay rates. Furthermore, we have demonstrated that this 
solution explains the new Belle data in the D^D'^n^ and k^k^J/ Xj/ decay modes and the n^n^J/ Xj/ BaBar data. The 
original Xci C^P) state acquires a significant DD* component and can be identified with the X (3940) . 
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